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Abstract 

In this paper we show that if A is a Poisson algebra equipped with a set of maps A^ ! ' : A — > 
A® N satisfying suitable conditions, then the images of the Casimir functions of A under 
the maps (that we call "loop coproducts") are in involution. Rational, trigonometric 
and elliptic Gaudin models can be recovered as particular cases of this result, and we 
show that the same happens for the integrable (or partially integrable) models that can be 
obtained through the so called coproduct method. 

On the other hand, this loop coproduct approach is potentially much more general, and 
could allow the generalization of the Gaudin algebras from the Lie-Poisson to the Poisson 
algebras context and, hopefully, the definition of new integrable models. 



1 Introduction 

The Gaudin model was defined in the '70 by Gaudin 12, 13J as an integrable quantum spin chain 
model. The model was later cast in the framework of linear r-matrix formulation by Sklyanin 
[19j and Jurco [15j . The linear r-matrix formulation associates classical integrable systems with 
loop algebras. In the case of the classical Gaudin models, given a simple Lie algebra g, one 
can define three different kinds of TV-body integrable sistems on the tensor product of N copies 
of the corresponding Lie-Poisson algebra ^"(fl*): rational, trigonometric and elliptic Gaudin 
models. To these integrable models correspond the rational, trigonometric and elliptic Gaudin 
algebras, respectively. If one considers the homogeneous rational Gaudin models, it turns out 
that some of its Hamiltonians are superintegrable, since they belong also to a set of involutivc 
functions obtained through the so called "coproduct method" [TJ [H H] and [5] for a review. 
The superintegrable Gaudin Hamiltonians are obtained when one applies the method using the 
standard coalgebra structure defined on a Lie-Poisson algebra thus establishing the relation 
between the coproduct method and the Gaudin models (see [3])- However, the second set of 
involutive functions generated through the coproduct method is complete only when g = sl(2). 
The problem of completing this set of functions has been solved when g = sl(N) in [10 using 
a bihamiltonian formulation. 

In the present paper we propose a generalisation of both the coproduct method and of the 
Gaudin algebras, that we call "loop coproduct method" . This generalisation works only in the 
classical case, but allows us to define a larger number of involutivc functions with respect to 
the standard coproduct method and, in particular, to recover all the missing integrals of the 



* Departamento de Fi'sica, Universidad de Burgos, E-09001 Burgos, Spain, email: fmusso@ubu.es 



1 



standard coproduct method when g is a simple Lie algebra. Moreover, rational, trigonometric 
and elliptic Gaudin models can be obtained as particular cases of this construction, so that 
it gives a unified theoretical framework for these algebraic approaches to classical integrable 
models. We also show how, through this generalisation, it is possible to treat in a natural way 
the rational Gaudin model with degeneracies and that, in this context, the coproduct method 
associated with the standard coalgebra structure of the Lie-Poisson algebras correspond to the 
case of maximal degeneracy. 

The plan of the paper is the following. In section [5] we briefly recall the coproduct method. 
In section [3] we recall the definition of the Gaudin models and of the Gaudin algebras. The 
loop-coproduct method is introduced in section [H where we also show that both the coproduct 
method and the Gaudin algebras are particular cases of this approach. In section[5]we show how 
rational Gaudin models with degeneracies can be obtained through the loop-coproduct method. 
Section O shows how applying this approach to Lie-Poisson algebras one can obtain a larger set 
of involutive functions with respect to the coproduct method. We derive in subsection 16.11 that 
when the Lie-Poisson algebra is simple these functions do indeed define a completely integrable 
system and in subsection 16.21 we give an explicit example for a non-simple algebra. Finally in 
section [7] we discuss the results we have presented and we comment on open problems. 

2 Coproduct method 

The coproduct method is a recipe to construct (classical and quantum) Hamiltonian integrable 
systems endowed with Poisson-coalgebra symmetry. In this section we tersely recall how it 
works in the classical case only, for detailed expositions, including different applications, we 
refer the reader to [U IH [7] . 

Let A be a unital, associative algebra and let us denote its identity map with id. If A is 
also endowed with a coproduct map 

A : A -> A ® A 

which is coassociative 

(A ® id) o A = (id ® A) o A 
and is an algebra morphism from A to A ® A: 

A(o6) = A(o)A(6), Va,beA, 

then (A, A) is called a coalgebra. If A is a Poisson algebra and A is also a Poisson morphism: 

A({o, b} A ) = {A(o), A(6)} A8A , Va, b e A, 

with respect to the standard poisson structure on A ® A: 

{a ® b, c <g) d} A ^ A = {a, c} A ® bd + ac ® {b, d} A a,b,c,d e A 

we shall say that (A, A) is a Poisson coalgebra. Let M be the dimension of A, {y a }a=i a se t 
of generators and let us suppose that r functionally independent Casimir functions Cj(y) = 
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Cj(y 1 , . . . ,y M ) are defined on A. Using the coproduct and the identity maps it is possible to 
construct the m-th coproduct maps: 



A (m) : A^A®A®--®A m = 2,...,N 

A (m) ;= ( A (m-1) g id j Q A (2) A^=A (1) 

that, in turn, are also both algebra and Poisson morphisms: 

A^iab) = A^(a) A™(b), Wa,beA, 



AW({a,i} i ) = (AW(a),AW(6)] Va,beA 



Let a be an element of A® A® ■ ■ ■ ® A and b be an element of A (8 A (8 • • ■ <S> A Then we 
define 

N—m N—n 



iV 



In the following all the Poisson brackets will be considered in A ® A ® • • • £g> A, so that the 
indication of the space in the Poisson bracket will be dropped. 

From the Poisson morphism property of the m— th coproduct it follows: 

{A( m )(C 4 ),A(")(y")} = if n>m. (2) 
From equation ([2]) it follows at once that: 

{A( ro )(C i ),AW(C i )} = 0, m,n = 2,...,N, i,j = l,...,r (3) 
{A(™) (Cj)jA W(y tt )}=0, m = 2,...,iV, i=l,...,r a=l,...,M (4) 
Hence, given a Poisson coalgebra A, the coproduct method allow us to find a large family of 

N 



functions in involution within the Poisson algebra A ® A Cg> • • • <g> A. 

In the following we will work mainly with the particular case when A is a Lie- Poisson algebra 
with commutation relations 

{y Q ,/} = C?V 

and the natural coalgebra structure 

A(y Q ) = y Q ® I + 1 <g> y Q = + . 
From now on, we will denote with {yf}, a = 1, . . . , M, i = 1, . . . , N the natural basis in 

N 

A<3 ■ ■ ■ (g) A: 

i-l N~i 

yf = 1 <g> • • • <g> 1 ®y a ® T" 

{//:'•///} ^ cs*y- 
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The m-th coproducts can be obtained from equation ([IJ) and in this basis they simply read 

m 

& (m) (y a ) = J2yi> ( 5 ) 
i=i 

then it is immediate to check that 

{a(*V),A«(/)} = C 7 ^A(*V) j>i. (6) 

Let us now go back to the general case of A being a Poisson algebra with r Casimirs 
function. It is natural to ask when the family of involutive functions obtained through the 
coproduct method will be large enough to define an integrable system. The dimension of the 
symplectic leaves of A is: 

, M -r 
d= —> 

N 

so the total number of degrees of freedom of an Hamiltonian system defined on A ® • • • A is 
N(M — r)/2. Equation defines a family of (N — l)r involutive functions, while equation 
(|1]) allow us to add to this family a maximal abelian subalgebra of A^ N \A) ~ A; that is at 
most (M — r)/2 further involutive functions. Hence, the coproduct method allow us to define 
an involutive family of at most 

M-r 
(jV_l) r +_— 

independent functions. Henceforth, a necessary condition for getting a complete integrable 
Hamiltonian system through the coproduct method is given by: 

M — r M — r 

(N - l)r + -— ^ > N — ^— => 2r > M - r, (7) 

and this condition is almost never satisfied for Poisson coalgebras. For example, if A is a Lie- 
Poisson coalgebra constructed out from a simple Lie algebra g, then the condition is satisfied 
only when g = sl(2). 

A natural question, when the coproduct method do not grant complete integrability, is if it 
is possible to enlarge in some way the family of integrals in order to define an integrable system. 
In section [fxll we show that, at least in the case of simple Lie-Poisson coalgebras, the answer 
to this question is affirmative. 



3 Gaudin models and Gaudin algebras 

In the modern formulation, TV-body classical Gaudin models are integrable models defined on 
the tensor product of N copies of a simple Lie-Poisson algebra ^(q*). Indeed, given a simple 
Lie algebra g and a solution r(A) G g ® of the Classical Yang-Baxter equation: 

[ri 3 (A),r 23 (/i)] + [ri2 (A - /i),r 13 ( A) +r 23 (M)] = 0, (8) 
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where 



ria(A) = r(A)®l, r 23 (A) = 1 ® r(A), n 3 (A) = (1 <8> n) ria(A) (1 8> II) 
and II is the permutation operator on g <g> g: 

n(a; ® y) = y <8> x, 

there exists a simple recipe (see, for example [17] or [20]) to construct a Lax matrix L(X) with 
entries in JF(g*)® ^ together with an r-matrix r'(A) satisfying: 



{L(A) 01,1® £(»} + [r'(A - /i), L(X) <8> 1 + 1 ® L(/i)] 



= 0. (9) 



Equation © ensures the commutativity of the spectral invariants of the Lax matrix L(X) that 
give rise to the Gaudin Hamiltonians. Drinfel'd and Belavin [5], proved that under suitable 
non-degeneracy conditions on r(A), the dependency on the spectral parameter can be only of 
three kinds: rational, trigonometric and elliptic. According to the chosen solution, equation 
([9]) will define the rational, trigonometric or elliptic Gaudin algebra (actually elliptic Gaudin 
algebras can be defined only for g = sl(n)). 

Let us make this construction explicit. Let p be a faithful representation of g, then the 
trace define a nondegenerate invariant bilinear form on p(g). Let {1"}*^, where M is the 
dimension of g, an orthonormal basis with respect to the trace, which means that in this basis 
the Cartan-Killing metric is diagonal 

Tr(X a X f3 ) = S at3 (10) 



a/3 _ 



fJ 

and the corresponding structure constants 

[X a ,X f3 ] = cfx^ 

are completely antisymmetric: 

In [5] it has been proven that the classical Yang-Baxter equation © always admits matrix 
solutions of the form (from now on we will write explicitly the sum in case Einstein convention 
could generate confusion): 

M 

r'(\) = Y J X a ®X a r(\) (11) 

a=l 

where /"(A) are suitable functions that can have an elliptic, trigonometric or rational depen- 
dence on the spectral parameter A. 

Let {y"}, a — 1, . . . , M be the coordinate functions on the Lie-Poisson algebra J r (Q*)'^ N cor- 
responding to the basis {X a }^ =1 . The "fundamental" Poisson brackets among the coordinates 
will then be given by: 

{y?,y?} = 5 ij Cfy7. 



1 For the sake of keeping the notation simple and make contact with the one used in [8] we decided to work in 
an orthonormal basis, hence considering g as a complex Lie algebra. The construction works "mutatis mutandis" 
also for the real forms of g and for an arbitrary choice of the basis (see 1201 ). 



5 



The Lax matrix of the N— body Gaudin model associated with the r-matrix is given by: 

M N 



a=l j=l 



(12) 



where tj are constant parameters subject to the constraint e; ^ tj, i 5^ j, the functions /"(A) 
are the same appearing in (jTTJ) and a is a constant matrix satisfying the constraint 



[a (g> 1 + 1 <g> cr,r(A)] = 0. 



(13) 



It can be proven (see for example [T7] or J5D]) that the Gaudin Lax matrix (|12p satisfies the 
r— matrix formulation ^ with the r— matrix pip . 

The Gaudin algebra generators are defined by 



JV 

y a (\) = Tr(L(\)X a ) = £/ Q (A - ej)yf + Tr(aX<*) 



(14) 



Let us evaluate the Poisson brackets between the Gaudin algebra generators ([T 

{y a (X), = Tr ({Z(A) ® 1, 1 ® £(M)}* a ® ^) = 

= -Tr ([r(A - m), L(X) ® 1 + 1 ® i(/i)]^ a ® X") = 



-Tr 



■£(i®xV(m) 



J2 / 7 (A - M)C £ ^Tr (((X e ® X V(A) + (Jf ® XVO*)) X" ® X* 



7<5 



= - ^ / 7 (A - ^)C7 tf (<P £ <PV (A) + S^S^'y 5 ^)) = 

7<5e 

= - E (A* " *0c£V(A) + r (a - ricfy'M) . 

s 

By using the fact that the structure constants are completely antisymmetric, finally we get: 
{y a (X), /(/,)} = -£ (/"(A - V(A) + /" (A - m)<^V (m)) - 

= E c f (/"( A - v)y s (p) - f(x - m)/(a)) . (15) 

s 

These are the defining relations of the Gaudin algebra associated with g and r' in the orthonor- 
mal basis (ITU1) . 
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3.1 The rational Gaudin model 



In the following we will deal mainly with the rational case: 

/ a (A) = i a = l,...,M. 
Accordingly, the Lax matrix of the rational Gaudin model associated with g is given by 

M N 



^) = E£f_r + - (16) 



a=l i=l 



In the rational case equation (fl3|) is always satisfied when a belongs to The rational 

Gaudin algebra generators are given by: 



N N 



y a (X) = Ti(L(X)X a ) = ' + Tr(aX a ) = £ + c a (17) 



, A - Ci f-f A - e, 

and the Gaudin algebra by the Poisson brackets: 

{y Q (A),/( M )} = -J—CfWW-ytM). (18) 
[i — X ' 

Let r = rankg and let us denote with C\,...,C r the fundamental Casimir functions in J-(q*). 
The rational Gaudin Hamiltonians are given by the spectral invariants of the Lax matrix (fH))) 
or, which is the same, by the residues of the Casimir functions evaluated on the rational Gaudin 
algebra generators (fTTjl : 

Tr {L(X) m ^) =C j (y\X),...,y M (X)) 

where rrij, j = 1, ...,r are the "exponents" [TT] of the Lie algebra g. When a has simple 
spectrum, the Gaudin Hamiltonians define a completely integrable system (provided that the 
parameters q are all distinct) 18J . If this is not the case, then there are missing integrals. 
These integrals can be recovered noticing that if r is a matrix belonging to the commutant of 
tr, i.e. [tr, r] = then the function: 

M N 



a—l i—1 

commutes with all the spectral invariants of (fH))) . By adding to the spectral invariants of (TIT))) 
a maximal abelian subalgebra of the Poisson algebra spanned by such functions 

M N 

{^ = ££yfTr(A-T) | [r,a)=Q}, (19) 

Q = l i=l 

complete integrability is recovered. 

When a — 0, the algebra (fTI?)) coincide with the whole algebra A^ N \j r (g*)), where A^' 
denotes the primitive iV-th coproduct map associated with the Lie-Poisson algebra F(g*) (see 



7 



equation [5]) . In such case, some of the rational Gaudin Hamiltonians turn out to be superinte- 
grable. Indeed, it can be shown that A( N \Ci), . . . ,A^ N ^(C r ), belong to the Gaudin Hamilto- 
nians [3 . Therefore, these Hamiltonians belong to two different sets of functions in involution, 
the one obtained from the spectral invariants of the Lax matrix (|16p and the one generated 
through the coproduct method (see equations [3] and |3J . 

4 Loop coproducts 

Let A be a generic Poisson algebra and {y a }^ =1 a set of generators for A with Poisson brackets: 

{y a ,y l3 } = F^(ff) y=(y\...,y M ). 
Let us denote with Y the set of all the generators on A® N : 

Y = {yf} a = l,...,M i = l,...,N 

i-l N-i 

yf = 1 ® ■ ■ ■ <8> 1 <8> y a ® 1 <g> • • • ® 1 

We state our main result: 

Theorem 1 Let us suppose that we have defined a set of m maps depending on a parameter X 

N 

A^ ] : A-» A®... A® A, k=l,...,m (20) 

and such that the following relations hold 

{Af( !/ a ) 1 AW(/)} = ^i,fc,A, /1 ,?)^(Af®) k>t (21) 
{Af (f),AW(/)} = 

= ^(fc,A, M ,f)F^(Ai fe) (y)) + ^(fc,A, M ,f)F^(A( fe )(y)) (22) 

for certain functions fS(i, k, A, fi, Y), gz(k, A, /Lt, Y), h"(k, A, fi, Y). 

If the map a(^ is defined on any smooth function of the generators f G A as: 

A«(/)( 2/ 1 ,...,^))= / (A«(y 1 ),...,A«(^)), 

then: 

{A«(C,),A«(/)} = fc>z (23) 
{A«(C,),AW(C0} = 0. (24) 

Proof: 
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Let us prove first the equation (|23p. We fix an arbitrary Casimir Cj. Since Cj is a Casimir 
function, for any /3 we must have: 

i c >y} = E S bf.v"} = E S*"*0> - ( 25 ) 

q=1 y a=l y 

Now we use this equation together with (j2"Tj) to prove equation ([25]) : 

M dC f A (l) (i/)^ 
{Af(Ci),Af(/)} = £ ^ {AiV), A«(/)} = 

= # (», fc, A, fi, Y) £ J (0 , / ^ 7 (Ai J) (y)) = fc > i 

a =i dA\'{y a ) 

From equation 

{A«(Cj),A<*V)} = fc>z 

it follows that 

{A»(C j ),AW(C I )} = z^fc 
so that only the case fc = i remains to be proven. In this case: 

{Al i) (C,-),A«(C i )} = 

= S « A (o, a / {a a V),a«(/)} 

By using equation (|22|) and (|25[) . we have: 
{A ( ^(C j ),A^(C l )} = 



aC, (A«(y)) ^ /m 9Ci ( A W(y)) , 

since the terms in parentheses vanish, which completes the proof. 



□ 



Hereafter we call the maps ([20)) as "loop coproducts" and the algebra ([2"T j) -(f22" |) as "general- 
ized Gaudin algebra" . In fact, it is straightforward to realize that the generators of the Gaudin 
algebras (see eq. (fT5jl ) are particular cases of (|2"T1) . (|2"2"j) provided that 

g^(k,X,^Y) = -5^f(X- fJ ,) 
h«{k,\ l x,Y) = ^r{X-ii)- 
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Also the coproduct construction is a particular case of Theorem [TJ since the m— th coproduct 
map satisfy (f2"Tj) and (f2"2")l for suitable choices of the arbitrary functions f@(i, k, A, fx, Y),gP (k, A, fx, Y) 
and hP{k,A,[i,Y). 

Let us show it. We recall, from section [5] that we have to consider the case of A being a 
Poisson coalgebra. Generically, the coproduct of a generator can be written as the finite sum 
(see W 

A^(y a ) = Y, G b(y)® H b(y) (26) 

b 

for some functions G£ , H". As shown in [2], thanks to the coassociativity property of the 
coproduct, the A:— th coproduct map can be written as 



A« = (a««A(^)oA< 2 ) i<k 
So that, putting together (|26p and ([27| . we have 

A {k) {y a ) =J2 G b( A(i) (y}) ® H^(A^(y)) i < k. 

b 

Let us now assume that i < k, by using (|28p we have: 

{A«V),A<*V)} =5] [{AW(y Q ),Gf(AW(y))}® H^A^M) 

b 



(27) 



(28) 



E 

6 L 

E 



<9AM(y7) 



-{AW(y a ),AW(yT)} Ufff(A^(^) 



which is just of the form (|21[) with 

*, A, fi, Y) = £ ® (AC*"*) W ) 

When k = i, from the homomorphism property of the coproduct, we simply have 

{A^(y a ),A^(y^)} = F^(A^(y)) 

which is of the form (12"2")) with (for example): 

g%(k,\,n,f)=6% h?(k,\,n,Y) = 0. 

Therefore, the loop coproduct construction can be interpreted as a simultaneous generalization 
of both the coproduct method and of the Gaudin algebra approach. 



10 



5 Rational Gaudin model with degeneracies 



In this section we will show how the case of the rational Gaudin model with degeneracies can 
be naturally formulated in the loop coproduct scheme. 

Let 2 be a simple Lie algebra as in section[3]and A — In the Subsection l3.1[ we have 

seen (see eq. (jT7]0 that the map 



given by: 



A A : A^ A® . . .A® A 



N 

A x (y a )^y a (\)=Y / T^ + c a (29) 



= 1 

with c Q , ei, . . . , ejv arbitrary constant parameters, defines the rational Gaudin algebra (| 1 8j) : 
{A x (y<*), A M (/)} = -^—Cf (A\(y y ) - A^)) 

/l — A 



Hence A\ satisfies the property (|22|) with 

1 



fi — A 

and the residues of the images of the Casimir functions under A\, A\(Ci), i = 1,. .. ,r define 
the rational g-Gaudin Hamiltonians: 

{A A (C i ),A M (C i )} = i,j = l,...,r (30) 

If the parameters e±, . . . , £n are all distinct: e; =/= €j, for i j and the constant matrix 

X a c° 



has simple spectrum, then equation pop provides enough functions in involution to get complete 
integrability. 

In the Subsection 13.11 we briefly discussed what happens when the spectrum of a is not 
simple. In this Section we will assume that a has simple spectrum and we will discuss the case 
when there are degeneracies in the parameters, i.e. e; = ej for some i ^ j. Also in this case, 
the spectral invariants of the Lax matrix (|16p are not enough to define a completely integrabile 
system. However, the Gaudin model with degeneracies is still integrable, and the expressions 
of the missing integrals was given in [TB] for the g — sl(2) case, and recently in [S] for the 
case of an arbitrary simple Lie-Poisson algebra. In the following we show that such missing 
integrals can be recovered through the loop coproduct formulation. Let us assume that we have 
k distinct values for the parameters ei, each with multiplicity di, i = 1, . . . , k. By renumbering 
the copies in ^4®^ and the e parameters, we can put 

£l = (2 = ■ ■ ■ = Cdi = Si 

£ cii+l = Erfi+2 = ' ' ' = Cdi+d, = S2 

. (31) 

EW-dfc + l = f-N-d k +2 = ■ ■ ■ = CjV = Sk 
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Let us introduce the k + 1 maps: 



Ai 2) (^)= E 



di a d 2 a N a 

ir ) ( „. )sW)=I :^ + v ^ + ... + v _i^ +c « 

where the r\i are iV distinct parameters such that r\i ^ rjj, i ^ j. These maps satisfy: 
{A«(y a ),A«(/)} = i,j = l,..., k 
{Af(f), A^V)} = A«(^) i = 1, . . . ,k 

{A^ +1 \y a ), A^ +1 \y f3 )} = -±^Cf (Aj^V) - A^V)) . 

By Theorem 1 we have that 

{Ai 4) (C0,A|/')(C m )} = i,j = l,...,A!+l, Z,m=l,...,r 

and by taking the residues of the functions A^(C;), complete integrability is recovered. If the e 
parameters are all distinct, then the residues of the function A^'(C;), i = 1, . . . , N give simply 
the Casimir functions of the i— th copy of A. 

Alternatively, we can complete the set of integrals in involution through a different set of 
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maps: 



3=1 

A^V^ E f + .- = 4 + 2,...,* 



Ar fe+1) (^)-W).EA-+ E A + E 



A — <5i ^— ' A — (^2 A — <5fc 

3=1 3=di+l j=JV-d fc +l 



Such maps satisfy: 



{a^V),a|/V)} = o 

ifdfe<« — k — 2 < dk+i , di < j — I — 2 < tf; + i ^ / 

{AlV),A«(/)} = icf A<V) 

if d fe < i - k - 2 < d k+ i . d k < j - k - 2 < d k+ i j > i 



l,...,N-k 



{AW(tf«),AW(^)} = (a£V) - A«(^)) 

f AW(tf a ), Af - fe+1 V)} = -^rCf A«(f) i = 1, . . . , TV - fc 
{Af- w V).Ar w V)} = (Af - fc+1 V) - Af "^V) 



(32) 



/' 

and they are also of the form required in Theorem 1. It follows that the images of the Casimir 
functions under this alternative set of maps give another family of involutive functions that, to 
the best of our knowledge, has not been considered in the literature. Since this two families 
have in common the Gaudin Hamiltonians 

A A (Ci), i = l,...,r 

it follows that the rational Gaudin model with degeneracies is superintegrable (see also P3]). 
The degree of superintegrability increase by increasing the degeneracy of the parameters e,, i = 
1,...,N. 
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6 Generalization of the coproduct method to higher rank 
Lie— Poisson algebras 



As we have remarked in section [5J the coproduct method is unable to provide complete in- 
tegrability when applied to a simple Lie-Poisson algebra with a rank higher than 1. In this 
section we will show how, using the loop coproduct formulation, it is possible to build up a 
larger set of involutive functions with respect to those obtained through the coproduct; in the 
next subsection we will show that, in the case of simple Lie algebras, these families are large 
enough to define completely integrable systems. 

Again, we consider the case when A is a Lie-Poisson algebra ^(g*). We recall that the 
primitive coproducts are given by (|6|): 

i 

A^(y a ) = J2vi i = 2 t ...,N 

j'=i 

Let us introduce the following loop coproducts 

Ai l V) = T^ + T^ i = 2,...,N (33) 

A A — C\ 

where e is a constant parameter (e ^ 0). Notice that these maps coincide with the first N — 1 
ones in (|3"2j) when k = 1 (that is when the parameters ej are all of them equal). The loop 
coproducts (|3"3")l satisfy the Poisson brackets: 

(A<V), A? >(/)} - + J«-, + Jt 

p L A A — e fx (X — 

= ~Cf ( + Y l ~)= V A £V) (34) 
fi ' \ A A — e J fi ' 

which are of the form (f!H) with f£(i, k, A, fx, Y) = 1/fx. 



, 7 . 

On the other hand, we have: 

{a£V),a£V)} 



A 



t V A/i ( M _ e )(A-e) 

- 1 gag f ^- 1) (y' y ) , a^H^) ^ j/ 



7 



/i — A 7 \ A A — e /i /i — 

^C^(a£V)-A<*V)) (35) 



which is of the form (22]) with g?(i, \, fj,,Y) = -h°(i, X, /jl,Y) = - A). Notice that ([35" 
coincide with the rational Gaudin algebra (fT8|) . From Theorem Q] it follows that: 

{A«(C J -)(^,AW(C / )(y)} = 
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Moreover, from (f!?5|) , it follows: 

{A«(C J -)(y),A«(y a )} = k>i Va (36) 

and, in particular 

{Al i) (C J -)(y),AW(y a )} = i = 2,...,JV Va. 

This result means that we can add to the involutive functions generated through the loop 
coproducts of the Casimirs A^ (Cj) a maximal abelian subalgebra of A^ N ^ (A) ~ A in the same 
way as in the standard coproduct case. 

Now let us show that the standard coproduct Hamiltonians are functionally dependent on 
the Hamiltonians obtained through the loop coproduct. The Casimir functions Cj(y) can be 
always taken to be homogeneous functions: 

C j (Xy\...,Xy M ) = f j (X)C j (y 1 ,...,y M ) (37) 

If we choose e = 0, then we have 

MX) - fj(x) MX) (2/)) " A {CAy)) 

So that we get back the standard coproduct Hamiltonians. Moreover, since 

{A ( £\c j )(y),A^(C J )(y)} = 

by virtue of (|37p . it follows that A 1 - 1 ' (Cj)(y) must be functionally dependent on the Hamiltonians 
obtained through the generating function A^(Cj)(y). 

6.1 Complete integrability in the case of simple Lie-Poisson algebras 

Let us show that through the loop coproducts (f3"3"| it is always possible to define a completely 
integrable system when A is a Lie-Poisson algebra ^(fl*) associated with a simple Lie algebra 
0- 

We recall from Section [5] that the number of degrees of freedom of an Hamiltonian system 

AT 

defined on A ® • • • <g> A is 

^ - (38) 

Let us denote with Cj the fundamental Casimirs of A, which are homogeneous polynomials 
in the generators of degree rrij + 1, where mj, j = 1, . . . , r are the exponents of the Lie algebra 

g pi]- Let us multiply A^(C^) by A"^ +1 , then 

A ™ J+ i A W (C 0)) = C U) U<.*-i)tf) + ^-yl . . . , A (i_1) (y M ) + ^-yf ) = 
= Af(C^)(y) 
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where 

Af (y°>) = A^- 1 ) (y a ) + \yf A = ^- (39) 
Let us expand A-^(C^) in terms of A 

Af (C«) = K^' j) + XK[ h3) + ■■■ + \ m i +1 K%f +1 (40) 

We notice that: 

Tnj + 1 



7(j) 

so that in the set of functions 



i = 2,...,N j = l,...,r 1 = 0,. 



) 3 



at most (N — 1) X)i=i m i are functionally independent among them and with respect to the 
Casimirs. From equation (|23[) it follows that we can add a maximal abelian subalgebra of 
A^ N \A) ~ A; that is at most (I — r)/2 further involutive functions. 

In conclusion, using the loop coproduct, we can define at most 

T 

(iV-l)£m j + (J-r)/2 
functions in involution. This set of functions will define an integrable system iff 

r 

2^ mj ->(Z-r). (41) 

3=1 

In the case of simple Lie algebras, formula (|4"Tj) means that twice the sum of all exponents 
has to be greater or equal to the dimension minus the rank. Let us denote with h the coxeter 
number associated with the simple Lie algebra, then the following two formulas holds: 

rh = l-r (42) 
mi + m r -i + i = h Vie(l, ...,r) (43) 

Summing (I43[) on all possible values of i and using equation (|4"2")) . we get: 

r r 

m, + m r ^i + i = 2 rrii = rh = I — r 



So that, for simple Lie algebras, the necessary condition for integrability (141|) is always verified. 
Moreover, it can be proven that the set of functions {K^'^} provide indeed (N — 1) x Y^j=i m j 
functionally independent involutive functions (see for example 18J), since the counting of in- 
dependent functions boils down to the counting of independent spectral invariants of the Lax 



1G 



matrix of the two-body homogeneous rational g-Gaudin model. In conclusion, in the case of 
simple Lie algebras, if we can define an integrable system on the corresponding Lie-Poisson 
algebra A, then, using the loop coproduct, we can define an integrable system on A® N . 

Finally, let us stress that an integrable system on a simple Lie-Poisson algebra can be always 
constructed through the spectral invariants of the following Lax matrix: 

L(X) = ^2x a y a + \a, 

a 

where a s p(g) must have simple spectrum. The iV-th coproduct of these Hamiltonians coincide 
with the residues of Ar^ f+1 ' (Cj), j = l,...,r, where A^ +1 ~* is the map in the last line of 
equation (|32p when k = 1. 



6.2 An example on a non-semisimple algebra 

We consider the two-photon Lie-Poisson coalgebra (he, A), that is a non-semisimple, rank two, 
six dimensional Lie-Poisson algebra spanned by the six generators {N, A + , A~ , B+ , B~, M} 
with Poisson brackets: 

{N,A+}=A+, {N,A~} = -A-, {A-,A+} = M, 

{N,B+} = 2B+, {N,B-} = -2B-, {B~, B+} = 2(M + 2N), 

{A+,B-} = -2A-, {A+,B+} = 0, {M, •} = 0, 

{A-,B+} = 2A+, {A-,B-} = 0, 

has two fundamental Casimir functions: the mass M and a third order Casimir given by 

C = M(N 2 - B+B- + MN) - A + A-(M + 2N) + B+(A-) 2 + (A+) 2 B-. 

The first step is to define an integrable system on /i 6 - Since is six dimensional with two 
Casimirs, we need two involutive functions. A (rather trivial) possible choice for involutive 
functions on h$ is: 

/i = B + + A~ + B~, h = (A+) 2 + (A-) 2 + (A-M). 
Notice that, from (|39|) . (140]) . At(C) will be a third order polynomial in A, the coefficients being: 
A : K (2) = C ® 1 

A 1 : K (2) = Mi(2NxN 2 - B+B- - B+B- + M X N 2 + M 2 N 1 ) + 
+ M 2 {N 2 + Mi Ni - B+ Bl ) - 

- A+A^(M 2 + 2N 2 ) - {A+Az + A+A^^Ah + 2N X ) + 
+ 2B+A-A- + B+(A-) 2 + (A+) 2 B- + 2A+A+B' 

X 2 : K { 2 2) = 1 <-► 2 

A 3 : K (2) = 1 ® C 

(2) (2) 

Kp has the same expression as K\ but with the indices 1 and 2 interchanged. By construction 
K^> and K^> are in involution among them and with A(/i) and A(/ 2 ). These four functions 
define an integrable system on hQ®h§. 
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We can define the following Poisson (iso-)morphism between (the symplectic leaves of) 
he ® he and the 8— dimensional symplectic algebra {qi,pj} — Sij i,j = 1, . . . , 4: 

D(A+ <g) 1) = MlPl + /i 2 P2 I>(1 <8 = M3P3 + M4P4 

® 1) = /iigi + /i 2 g2 -0(1 ® A - ) = ^3(73 + /i 4 g 4 



£>(M ® 1) = Mi + I>(1 ® Af) = A*l + A»4 
o\+pl+- ( ~ V2 D(l®B+)=pl+pl + - 5 

£>(£" <g> 1) = g? + gl D(l®B-)=gf + g* 



D(JV ® 1) = p l5l + p 2 g 2 - D(l ® TV) = p 3 q 3 + p m - 

With this choice the functions 

{D(A(h)),D(A(f 2 )), D(K{ 2) ), D{K { 2 2) )} 
define an integrable system on the 8— dimensional symplectic algebra. In particular, we have: 

^ (M192-M291) (^394-^453) 

D(A(f 2 )) = (l>^) + (^^) 

(2) (21 

The explicit expressions for D(K[ ') and D(K 2 ') are really cumbersome, so that we do not 
write therr0. 



7 Conclusions 

In this paper we have presented a generalization of both the Gaudin algebras and the coproduct 
method that we called loop-coproduct method. We showed that this method provides a unified 
theoretical framework for rational, trigonometric and elliptic Gaudin models, and also for the 
integrable models coming from the coproduct method. Moreover we showed how, in the case 
of Lie-Poisson algebras, the loop-coproduct method makes it possible to define larger sets of 
involutive functions with respect to the coproduct method and that this larger sets are indeed 
complete when the Lie-Poisson algebra is simple. In this section we would like to mention 
some open questions that, in our opinion, deserve further investigation. One of the questions 
is whether this approach can be extended also to quantum systems, that is how one should 
reformulate theorem[T]whcn A is a noncommutative algebra (with respect to the multiplication) . 
The quantization of the rational Gaudin Hamiltonians, found in [21j . could be a starting point 
to investigate this problem. 

2 In [7] a degenerate symplectic realization has been considered. In this case, to obtain an integrable system, 
one can indeed consider only the standard coproduct of the Casimirs but, at the same time, still has to begin 
with an integrable system defined on the abstract algebra hg. 



18 



The only examples of integrable systems obtained from theorem [T] that we discussed in 
this paper are related to Lie-Poisson algebras. On the other hand the theorem is much more 
general, since can be applied to arbitrary Poisson algebras. Indeed, as the coproduct method 
is a particular case of theorem [1] the integrable systems coming from the coproduct method 
applied to q-Poisson algebras (see [5]) also fall into this framework. However, these are the 
only examples that we know of integrable systems that can be obtained by applying the loop- 
coproduct method to a (non-linear) Poisson algebra and in all these examples there is no 
dependence on the "spectral" parameters A and \x. Equation (|22[) can be seen as a defining 
equation for A^ . Indeed, it can be rewritten in the form 

v aAl%°)aA<, fc V) , 

= g?(k, A, fi, Y)F a -<(A[ k \y)) + h»(k, A, fi, Y)F^(A^(y)) (44) 

This is a non-linear PDE with a functional part in the spectral parameters A, fi. Clearly, finding 
a solution (if it exists) is a formidable task. On the other hand, in equation (|4~4"|) , we have the 
freedom of choosing the Poisson algebra (that is the functions F a P) and the arbitrary functions 
gP and h". A clever choice of such functions could make the solution of equation (|4*4")) affordable. 
In our opinion, it would be very interesting to find other maps defined on a Poisson algebra, 
depending on the spectral parameters and satisfying the hypotheses of the theorem [TJ Indeed, 
such maps would define a natural generalisation of the Gaudin algebras from the Lie-Poisson to 
the Poisson algebra context and could, hopefully, provide examples of new integrable models. 
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